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NATIONAL  ADVISORY  COMMITTEE  FOR  AERONAUTICS 


TECHNICAL  NOTE  NO-  1852 

REINFORCED  CIRCULAR  CUTOUTS  IN  PLANE  SHEETS 
By  H-  Reissner  and  M-  Morduchow 


SUMMARY 


The  problem  treated  here  is  to  design  the  reinforcement  of  a  cutout 
in  a  plane  sheet  in  such  a  way  that  it  is  as  nearly  as  possible  equivalent 
to  the  part  of  the  structure  which  has  been  cut  out.  A  perfect  equivalence 
would  mean  that  the  stresses  and  displacements  of  the  structure  remain  the 
same  as  those  which  would  have  appeared  without  the  cutout. 

General  formulas  are  developed  for  the  circumferential  distribution  of 
the  cross-sectional  moment  of  inertia  Ir  and  of  the  area  Ar  of  a 

circular  reinforcement  required  for  perfect  equivalence.  These  formulas 
are  then  applied  to  some  cases  of  external  edge  tractions;  Hydrostatic 
stress,  pure  shear,  uniaxial  tension,  and  pure  bending.  It  is  found  that 
in  the  first  two  cases,  the  required  cross  sections  are  physically  possible 
(i.e.,  Ir  and  Ar  come  out  positive),  although  the  required  moment  of 
inertia  is  in  some  cases  found  to  be  quite  high  in  comparison  with  the 
required  area. 

In  the  cases  of  uniaxial  tension  and  pure  bending,  it  is  shown  that 
constraint  stresses,  that  is,  additional  stresses  in  the  sheet  due  to  the 
reinforced  cutout,  are  practically  unavoidable.  Simple  formulas  are 
developed  for  calculating  these  Constraint M  stresses  for  any  given  (constant) 
cross-sectional  characteristics  of  the  reinforcement  ring.  These  formulas 
are  derived  on  the  basis  of  the  assumption  that  the  constraints  diminish 
sufficiently  rapidly  with  radial  distance  from  the  cutout  so  as  to  have 
little  effect  at  the  external  edges  of  the  sheet. 

To  check  the  influence  which  actual  boundary  conditions  might  have  on 
the  practical  validity  of  these  formulas,  a  test  was  made  on  a  plane  sheet 
with  a  reinforced  circular  cutout  subjected  to  a  tensile  load  causing 
constant  displacements  at  the  loaded  edges.  It  was  found  that  the  values 
of  the  strains  calculated  from  the  exact  formulas  developed  in  this  report 
for  an  infinite  sheet  were  fairly  similar  to  the  values  of  the  measured 
strains  in  the  specimen,  except  along  the  loaded  edges,  where  the  actual 
strain  decreased  more  rapidly  toward  the  center  of  the  sheet  than  the 
calculated  strains.  This  discrepancy  must  be  due,  at  least  in  part,  to 
the  actual  condition  of  constant  displacements  at  the  loaded  edges  of 
finite  length  instead  of,  as  in  the  analytical  formulas,  constant  stress 
at  the  remote  loaded  edges  of  infinite  length. 
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INTRODUCTION 


This  investigation  is  concerned  with  the  problem  of  the  reinforcement 
of  cutouts  in  plane  sheets.  Such  cutouts  ma; y  serve  either  to  make  a 
structure  lighter  or  to  provide  space  for  personnel  or  accessories.  The 
problem  of  designing  a  reinforcement  ring  around  the  cutout  so  that  it  will 
be  elastically  equivalent  to  the  portion  of  the  sheet  removed  is  treated 
here  by  an  inverse  method.  The  treatment  is  given  first  in  general  terms 
and  then  applied  in  detail  for  several  types  of  edge  loading. 

A  test  was  made  on  a  plane  sheet  with  a  reinforced  cutout  to  check 
the  effect  of  seme  of  the  simplifying  assumptions  on  the  formulas  developed. 

This  investigation  was  carried  out  at  the  Polytechnic  Institute  of 
Brooklyn  under  the  sponsorship  and  with  the  financial  assistance  of  the 
National  Advisory  Committee  for  Aeronautics.  Grateful  acknowledgement  is 
hereby  also  made  to  Professor  N.  J.  Hoff  and  to  Dr-  B«  Boley  for  giving 
the  authors  the  benefit  of  the  experience  of  the  Laboratory  for  Aircraft 
Structures  of  the  Polytechnic  Institute  of  Brooklyn. 


SYMBOLS 


A© 

Ar 

a 

al 

b 

cmn 

d 


e 


Gr 
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effective  area  of  reinforcement  ring  with  rivet  holes 

cross-sectional  area  of  cutout  reinforcement  ring 

radius  of  center  line  of  ring 

radius  of  outer  circumference  of  ring 

half  width  of  sheet 

arbitrary  constants  in  stress  function 
radial  width  of  ring 

distance  from  edge  to  neutral  fiber  of  reinforcement  ring 

modulus  of  elasticity  of  ring  material 

modulus  of  elasticity  of  sheet  material 

bending  moment  in  cross  section  of  ring 

height  of  web 
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moment  of  inertia;  with  subscripts  e,  F,  and  Vy 

the  effective  moment  of  inertia,  moment  of  inertia  of 
the  flanges,  and  moment  of  inertia  of  the  web, 
respectively 


cross-sectional  moment  of  inertia  of  ring 


dimensionless  parameter  of  moment  of  inertia 


half  length  of  sheet 

abbreviations  for  products  of  and  powers  of  a^ 

bending  moment  about  z-axis 

radial  shear  force  on  cross  section  of  reinforcement  ring 
order  of  terms  of  stress  function 


load 

coefficients  in  stress  function  (cf.  equation  (3a)) 


radial  distance  from  center  of  cutout 

radial  distance  from  center  of  cutout  to  neutral  fiber 
thickness  of  sheet 
normal  stress  resultant 

displacements  in  x-  and  y-directions,  respectively 

same  displacement  system  in  radial  and  transverse 
directions,  respectively 

Cartesian  coordinates,  measured  from  center  of  cutout 

tangential  and  radial  loads  on  ring  per  unit  of 
circumference 

(E  Ax* 

shear  strain 

dimensionless  parameter  of  width  of  ring  (d/2a) 
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strain  in  x-  and  y-directions,  respectively 
radial  and  transverse  strain,  respectively 
change  of  curvature  of  ring 
Poisson’s  ratio 
radius  of  curvature 

normal  and  shear  stresses  in  sheet  (Cartesian  coordinates) 

radial,  transverse,  and  shear  stresses  in  sheet  (polar 
coordinates) 

original  shear  stress 

angular  polar  coordinate  (fig*  l) 

stress  (Airy)  function 


Suh - subs cr i pt s : 
o  original 

1  constraint 

T  total 

GENERAL  EQUATIONS  AND  GENERAL  SOLUTIONS 


In  order  to  analyze  the  effect  of  a  circular  cutout  in  a  plane  sheet, 
the  stresses,  displacements,  and  strains  in  such  a  sheet  may  he  represented 
in  a  polar-coordinate  system  (fig-  1)*  The  original  stresses  and  displace¬ 
ments,  before  being  disturbed  by  the  cutout,  will  then  appear  in  the  form 

\ 


(1) 
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The  cutout  edges  are  assumed,  to  carry  reinforcement  strips 
is  the  purpose  of  the  analysis  of  this  paper  to  choose  the 
of  such  an  edge  reinforcement  so  that  the  increments  or  , 


or  angles.  It 
elastic  properties 


and 


% 


to  the  original  stresses  (the  original  stresses  heing  given  by  the  sheet 
edge  loading)  hecome  zero  or  as  small  as  possible. 


Stresses 

The  stress  in  a  plane  sheet  can  he  expressed  by  a  stress  function 
\|r(r,  cp)  such  that  *'') 


o  =  r’^ty  4.  \|r * 

%  -  *n 
Tr9'-(r"1+)' 
where  the  primes  and  dots  are  defined  by 


(2) 


1  ^  . 

Sr 

S  = 


The  function 


^ ,  moreover,  must  be  such  that  it  satisfies  the  equation 


/V\t[f  =  0 


where 


Ail r 


S2 1  S2 1 

dx2  dy2 


_  1  d  /  d\|r\  ,  1  d2i|r 
“  r  dr  \  dr/  r^  dq£ 


(3) 


The  complete  solution  of  equation  (3)  for  the  case  considered  here, 
in  which  the  edge  stresses  taken  around  the  circumference  of  the  circular 
cutout  have  zero  resultant,  may  be  written  in  the  form: 
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t  Pn  cos  nqp  +  sin  ncp 

n=0  n=l 


(3a) 


where 


\ 


P0  “  °10  *  °20r2  +  c30  lc«9  r  +  So1-2  ±0%  r 


P1  •  cllr  *  c21r3  +  +  °lar  lo®e  r 


Pni2  =  °mrn  +  c2nr_n  +  +  c'-r 


n+2  -(n-2) 

3n  C4nr 


(3b) 


J 


The  expressions  for  the  functions  Q  are  the  same,  except  with 
possibly  different  numerical  values  of  the  constants  c- 

From  equations  (2)  and  (3a)  it  follows  that 

°r  ”  ^  (-n2r~2pn  +  r“1Pn1^  cos  ncp 

n=0  '  7 

-Yp" 

qp  /  rn 


=  /  rn  cos  ncp 
n-0 


w 


rep 


=  ^r”-1-?^)  n  sin  ncp 


n-0 


The  terms  in  Q  may  be  obtained  from  equations  (4)  by  exchanging  cosine 
for  sine  and  vice  versa  and  by  making  the  right-hand  side  minus  in  the 
last  equation  for  Tr(p* 


Displacements 


The  radial  and  transverse  strains 
from  the  stresses  by  Hooke fs  law: 


€r  and  e  in  the  sheet  follow 
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®aer  “  ■Esur  -  ar  ~  VcJcp 


E  e 
s  cp 


\ 


> 


) 


(5) 


where  Eg  is  the  modulus  of  elasticity  of  the  sheet  material. 

The  stress-strain  relations  of  the  reinforcement  ring  can  he  expressed 
in  the  form  ( see  appendix  A) : 


era  = 


Ta 

EpAj. 


+  u.  — 


2a 


(6a) 


Ga 

ErJr 


=  Up  +  Up 


(6b) 


where,  for  the  ring,  €r  and  *r  are,  respectively,  the  extension  strain 

and  curvature,  T  is  the  normal  stress  resultant,  G  the  bending  moment 
(at  any  cp) ,  d  the  width  of  the  ring,  Ap  and  Ir  the  cross-sectional 

area  and  moment  of  inertia  of  the  ring  (which  may  be  functions  of  cp) , 
and  U|_  and  Up  the  circumferential  and  radial  displacements  of  the  ring 

at  its  outer  circumference  (r  =  ap) . 


Since  the  displacements  of  the  ring  must  be  the  same  as  those  of  the 
sheet  where  the  ring  is  joined  to  the  sheet,  it  follows  that  the  displace¬ 
ments  in ’equations  (6a)  and  (6b)  can  be  obtained  by  putting  r  =  ap  in 
equations  (5)  for  the  displacements  in  the  sheet. 


The  stress  and  moment  resultants  T  and  G  in  the  ring  can  be 
derived  from  the  equations  of  equilibrium  of  the  ring  with  the  use  of  the 
fact  that  the  unit  loads  acting  on  the  reinforcement  ring  are  due  to  the 
radial  and  shear  stresses  ar  and  t  in  the  sheet  along  the  circumference 


of  the  cutout  ring. 


rep 


Cross-Sectional  Properties  in  Terms  of  the  Airy  Function 

Substituting  the  expressions,  as  obtained  in  the  preceding  discussion 
(see  appendix  B) ,  for  the  displacements  and  the  stress  resultants  into 
equations  (6a)  and  (6b)  and  solving  these  equations  for  Ap  and  Ir  yield 


the  following  expressions  for  the  cross-sectional  area  and  moment  of  inertia  of  the  cutout 
reinforcement  ring: 
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and  Qn  into  equations  (7a)  and  (71))  gives  the  distribution  of  cross-sectional  cutout-ring 

area  and  moment  of  inertia  required  to  give  zero  constraints,  that  is,  required  to  replace 
with  el.astic  equivalence  the  portion  of  the  sheet  cut  out*  If  the  values  of  Ar  and  Ir  thus 
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obtained  are  positive  for  any  cp,  then  these  are  the  correct  and  only- 
values  for  zero  constraints.  If,  on  the  other  hand,  the  resulting  values 
of  Ar  and  Ir  are  negative  or  infinite  for  certain  values  of  cp,  then 

an  elastically  equivalent  reinforcement  cannot  be  obtained.  This  procedure 
is  considered  in  more  detail  in  the  following  discussion. 


In  equations  (7a)  and  (7b),  the  Pfs  (or  Q's)  in  the  general  case 
represent  the  total  stresses,  that  is,  original  plus  "constraint"  stresses, 
in  the  sheet.  In  case  the  reinforcement  ring  around  the  cutout  is  equi valent 
to  the  cutout  portion  of  the  sheet,  then  the  constraint  stresses  will  be 
zero  and  only  the  original  stress  distribution  will  exist  in  the  sheet. 

Hence  In  that  case,  and  in  that  case  only,  the  P's  will  represent  the 
original  stresses.  The  original  stresses  in  a  sheet  due  to  any  loading  of 
the  external  edges  must  be  finite  at  all  finite  values  of  r  including  the 
center  (r  =  0) ;  hence  they  can  be  represented  in  general  by  an  Airy  func¬ 
tion  of  the  form  given  by  equation  (3a)  in  which  the  coefficients  Pn 
(or  On)  have  the  form  (cf.  equations  (3b)): 


P0  =  c20r‘ 


P1  °  °21r3 


P  =  c  Tn  +  c  rn+^ 


(3c) 


n~2  In 


3n 


The  expressions  (3c)  can,  if  dSsired,  be  put  into  equations  (7a)  and  (7b)  to 
give  slightly  more  explicit  expressions  for  Ar  and  Ip  (as  functions 

of  cp)  in  terms  of  the  constants  c  which  depend  on  the  given  edge-loading 
conditions. 


SPECIAL  CASES 

Cases  of  Zero  Constraint  Stresses 


The  problem  of  the  choice  of  an  exactly  equivalent  reinforcement  ring 
around  a  circular  cutout  in  a  plane  sheet  Is,  in  principle,  completely 
solved  in  GENERAL  EQUATIONS  AND  GENERAL  SOLUTIONS.  It  remains  to  show  some 
of  the  practical  implications  of  this  solution.  For  this  purpose,  the 
special  case  will  first  be  considered  in  which  the  edge  loading  on  the 
sheet  is  such  that  the  original  stress  distribution  can  be  represented  by 
an  Airy  function  with  only  a  single  trigonometric  term,  that  is. 
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i|r  =  P„  cos  ncp 


t  =  On  sin  n(P 


where  Pn  and  are  given  "by  expressions  (3c) •  This  class  of  cases 

will  lead  to  cross  sections  constant  along  the  circumference  of  the  ring. 


For  n  =  0,  equations  (7a)  and  (7h)  (c2  =  °3  =  ^  ^ecause 

radial  symmetr^  lead  to  the  result: 


of  the 


^r  = 


at  *_e 
1  -  V  E 


Tr  2(1  -  V) 


2a  ate- 


The  fact  that  the  constant  Cj_  is  arbitrary  shows  that  in  this  case  the 

moment  of  inertia  I  is  arbitrary;  but  the  cross-sectional  area  must, 
for  zero  constraints,  have  the  value  given  by  equation  (8a) . 

For  n  =  1  there  results  (assuming  symmetry  with  respect  to  the 
y-axis  and  setting,  therefore,  C2  =  0) : 


Ar  = 


k2i  [6  -  2V  -  ^  (1  -  3V)] 


2kon  0-  +  -  k  cos  9  +  — 

Eq  21V  2a/  3  at 

I„  =  a3t-2  ± - — - 


vhere,  for  abbreviation. 


-  c3 
k3  t 


k21  =  c2ia 
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In  order  that 


Ir  he  finite,  it  is  necessary  that 


Whence , 


(9) 


Again,  as  for  the  case  n  =  0,  Ir  may  be  chosen  as  desired,  but  the 
area  Aj,  required  for  zero  constraints  is  fixed  by  equation  (9). 

For  n  >  2,  by  setting  c-^  =  c2  =  c^  =  0  to  obtain  a  constant  cross 
section,  equations  (7a)  and  (Tb)  lead  to: 


Es 

Ar  =  at^ 


Kh  +  n  +  2 


K 


n(n-i.  1, 


n  (1+  V) 


Ir  =  aJt 


3*  ® 


E, 


1  k  -  x)k  -  j +  [(n  +  x)t  -  j 

n2  -  1  -n(l  +  v)k  +  [2  -  n  -  v  (2  +  nj] 


(n  +  l)  jjn  -  2)v  +  n+ 2J  +  ^  j^(n  +  2)v+n-2 

(10a) 


(10b) 


where 


_  cln 

K  =  — — 


c  a 
3n  1 


and  K  is  prescribed  by  the  original  stress  distribution.  Here  both  the 
<  moment  of  inertia  and  the  area  are  fixed  if  zero  constraints  are  to  be 

achieved. 
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Two  technically  important  cases  of  exactly  equivalent  reinforcements 
having  constant  cross  sections  are  included  in  equations  (8a) ,  (8h) ,  (10a) , 
and  (10b)  and  are  those  of  hydrostatic  stress  and  of  pure  shear.  This  can  * 

he  readily  seen  in  the  following  discussion. 

Homogeneous  hydrostatic  stress  ar  =  a^,  TrCp  =  0.-  In  this  case  the 

Airy  function,  as  may  he  verified  by*  equations  (2) ,  is  given  by  equations  (3b) 
as 


*  =  P0  =  c20p2 
°r  =  acp  =  2c20 

From  equations  (7a)  and  (7b)  it  follows,  as  in  equations  (8a)  and  (8b),  that 

at  >fl 
r  ’  1  -  V  Er 

Ir  =  Arbitrary  value 


(11a) 


(Hb) 


These  equations  show  that  the  equivalence  of  the  reinforcement  is  assured  for 
the  definite  value  of  the  cross  section  Ap  given  by  equation  (lla)  but  for 

an  arbitrary  value  of  the  moment  of  inertia  Ir  of  the  ring.  The  value 

of  Ir  will  then,  of  course,  be  chosen  as  small  as  is  compatible  with*  Ar 
and  with  buckling  considerations. 

Homogeneous  shear  distribution. -  This  (original)  stress  distribution 
can  be  expressed  by  TXy  =  To  in  Cartesian  coordinates.  In  polar 

coordinates, 

(12) 

J 


oT  =  -tq  sin  2cp 

c_  =  r  sin  2 cp 
cp  o  Y 

T  =  ““T-  COS  2cp 

rep  0 


Ir' 


/ 
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The  corresponding  Airy  function,  from  equations  (2),  is  seen  to  he 


*  = 


T 

O 


~  sin  2<p 


(12a) 


signifying  that  (cf.  equations  (3c)) 


Qq  —  c^gr 


where 


C12  "  To/2 


c32  “  0 


Using  equations  (7a)  and  (7b),  but  with  sine  terms  instead  of  cosine  terms 
and  with  cp  =  C2  =  Cg  =  0  to  obtain  a  constant  cross  section,  gives 


Ar  - 


at 


(1  +  v)(l  + 


4dN 

2a, 


a3t 


E, 


6(1  +  v)  ®, 


(13) 


(14) 


States  of  Stress  with  Imperfectly  Equivalent  Reinforcement  Rings 

In  the  previous  special  cases  of  one  trigonometric  term  for  the  Airy 
function  corresponding  to  the  original  stress  distribution  in  the  sheet,  it 
was  seen  that  a  reinforcement  around  a  cutout  could,  theoretically,  be 
designed  so  as  to  produce  no  additional  stresses  due  to  the  cutout.  The 
required  cross  sections  were,  in  fact,  constant  along  the  circumference. 

For  most  other  cases  of  external  edge  loading,  however,  it  will  be  found 
that  equations  (7a)  and  (7b),  intended  to  give  full  equivalence,  will  lead 
to  physically  impossible  (i.e.,  negative  or  infinite)  values  of  A^  and  Ir 

for  at  least  some  values  of  cp  and  that  zero  constraints  are  therefore 
impossible.  Moreover,  even  in  cases  where  zero  constraints  can  be  theore¬ 
tically  achieved,  it  may  be  found  (as  in  the  case  of  pure  shear,  discussed 
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further  in  the  following  section)  that  the  required  geometric  properties 
of  the  cross  sections  are  in  practice  not  realizable.  In  such  cases  it  is 
of  interest  to  determine  how  small  the  constraint  stresses  can  be  made  for 
(preferably)  constant  cross  sections  of  a  reinforcement  ring. 

The  constraint  stresses  can  be  calculated  without  difficulty  by 
substituting  some  given  constant  values  of  A^.  and  Ir  into  equations  (7a) 
and  (7b) ,  using  expressions  (3b)  for  the  Pn  (or  Qn) ,  including  both  the 
original  (given)  and  constraint  (unknown)  terms,  and  then  determining  the 
unknown  constants  c2n  and  c^n  (corresponding  to  the  constraint  stresses) 
so  that  equations  (7a)  and  (7b)  are  identically  satisfied  for  any  angle  cp. 

The  constraint  stresses  will  be  calculated  here  under  the  assumption  that 
they  vanish  at  infinity.  A  number  of  constants  then  disappear  in  the 
expressions  (3b).  This  condition  is  selected  to  give  practically  negligible 
constraint  stresses  near  the  edges  of  the  sheet*  Because  of  this  condition 
the  constraint  stresses  in  the  sheet  as  derived  in  this  report  must  be 
considered  as  approximate  (they  would  be  exact  for  an  infinite  sheet),  but 
the  approximation  will  be  good  if,  as  is  commonly  the  case,  the  constraint 
stresses  diminish  sufficiently  rapidly  with  distance  from  the  cutout  so  as 
to  have  a  very  small  or  negligible  magnitude  at  the  edges  of  the  sheet. 

Three  special  cases  of  edge  loading  are  now  considered  in  detail. 

Pure  shear.  -  It  will  be  found,  upon  closer  examination  of  expressions  (13) 
and  (l4),  that  the  moment  of  inertia  of  the  ring  required  for  an 
exactly  equi valent  reinforcement  in  the  case  of  pure  shear  Is  very  high 
if  its  cross-sectional  area  be  not  higher  than  that  required  by  equation  (13)« 
In  actual  design,  therefore,  the  cross-sectional  properties  given  by 
equations  (13)  and  (14)  cannot  ordinarily  be  realized,  and  constraint 
stresses  must  be  allowed.  In  the  following  paragraphs  the  method  of  cal¬ 
culating  these  stresses  is  given  in  detail. 

For  pure  shear  it  suffices,  except  for  the  addition  of  terms  in  r, 
to  employ  only  the  trigonometric  term  in  the  series  for  the  Airy  stress 
function  t  which  was  used  in  the  case  of  no  constraints,  namely: 


i|r  =  Q2  sin  2<p 


where  (cf.  equations  (3b)) 


-  cl2r 


2  + 


c22r 


-2  +  c 


42 


and,  as  before , 

_T_c 

C12  2 


(15) 


V 
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The  constraints  are  given  "by  the  terms  in  C22  and  which  are  the 

unknowns.  In  expressions  ( 3h) ,  the  constant  c^g  must  he  made  zero  because 

otherwise  it  would  give  constraint  stresses  which  do  not  decrease  with 
distance  from  the  cutout.  Putting  equation  (15)  into  equations  (7a)  and  (7b) 
and  introducing  for  abbreviation 


k22  =  022*1  > 

fcMaVai'2  , 

the  expressions  for  A^.  and  1^.  become: 


(15a) 


E, 


at 


t  o  ""  21^22 


A  —  — 
r  Er  1  +  V 


T  (1  +  kb)  +  (6  -  85 )k 


22 


/  165 
\1  +  V 


(16a) 


3v 


1  +  v. 


k42 


E, 


a3t 


to(1  -  5)  +  (2  +  6s)k22  +  (2  +  25)k^2 


Ir  Er  6(1  +  V) 


2k22 


_i_k 

1  +  V 


(16b) 


Let 


^r  at 

Es  ■r  "  2(1  +  v)a 

Ep  _  a^t 
^Ir  =  24(1  +  V)1 


Then  solving  equations  (l6a)  and  (l6h)  for  k22  ^42  with  v  =0-3  gives 
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of  equations  (17) ,  and  the  total  stresses  in  the  sheet  can  then  he  calculated  hy  determining  the 
constraint  stresses  as  given  hy  equations  (l8)  and  adding  these  to  the  original  stresses, 
equations  (12). 
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The  following  table  gives  the  values  of  k22  and  k^2  for  several 

reasonably  practical  sets  of  values  of  thickness,  area,  and  the  stiffness 
ratios  5,  a,  and  i.  The  entry  values  for  5  =  0  must  be  understood  to 
be  an  approximation  for  a  very  small  width  of  the  ring.  The  approximation 
i  =  0  for  very  small  values  of  5  introduces  only  negligible  errors. 


PURE  SHEAR 
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It  will  be  seen  fr cm  this  table  and  equations  (l8)  and  (12)  that  in 
actual  design  it  is  difficult  to  avoid  entirely  some  stress  concentration 

at  the  cutout  ( —  =  lV  Without  any  reinforcement  (i.e.,  for  &=i=a= 

V*1  / 

the  stress  concentration  at  the  cutout  is  extremely  high,  for 


where  a  is  the  total  tangential  normal  stress  with  cutout  and  0  is 
CP*p 

the  original  stress  (without  cutout) .  It  will  he  found  that  reinforcements 
with  practical  cross  sections,  such  as  those  in  the  foregoing  table,  will 
relieve  this  normal-stress  concentration,  although  they  will  introduce  a 
smaller  shear-stress  concentration  J’rom  the  point  of  view  of  minimum 

ratios  of  total  stresses  to  local  original  stresses,  the  most  suitable 
cross  section  given  in  the  foregoing  table  is  5  =  0.1,  i  =  0.1,  a  =  1. 

The  stresses  for  this  section  are  given  in  the  following  table.  (Also, 
see  table  I. ) 


HIRE  SHEAR 


at  _  aot  .  ..  ^  ^  \ 

5  ■  °'1;  '  H T77)'  Ir  ‘  240(1  -  V  ‘  °-3J 


r/a  1 

%r° 

Q 

o 

Trq >l/T° 

°rT/0(Po 

crqjr/acpo 

Tr^lTr% 

1.0 

l.l (-6 

0.303 

0.577 

-0.46 

1-303 

1-577 

1.25 

1.00 

.124 

.438 

.00 

1.124 

1.438 

Uniform  single  tension  or  compression.  -  In  Cartesian  coordinates,  this 
state  of  original  stress  is  given  by 


T 

^O 


=  0 


where  o  is  a  constant, 
o 
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In  polar  coordinates,  the  same  state  of  stress  is: 


orQ  =  -]p(2  -  2  cos  2cp) 


=  -jj-(2  +  2  cos  2cp) 


T  =  -ip  2  sin  2cp 
r<Po  4 


> 


(19) 


y 


The  Airy  function  following  from  equations  (2)  and  (19)  has  the  form. 


t  =  c2Qr2  +  c12r2  cos  2cp 
Thus, 

=  -^r2(i  +  cos  2<p)  (20) 

It  is  easy  to  show  that  in  this  case,  a  perfectly  equivalent  reinforce¬ 
ment  ring  is  not  realizable  (see  appendix  C).  The  Airy  function  for  the 
final  state  of  stress  in  the  sheet  will  therefore  contain  terms  representing 
stresses  due  to  the  reinforced  cutout,  that  is,  constraint  terms,  in  addi¬ 
tion  to  those  corresponding  to  the  original  state  of  stress.  Beferring  to 
equation  (20) ,  it  is  seen  that  the  Airy  function  for  the  original  (super¬ 
script  o)  state  of  stress  is  of  the  form 


Here, 


c20  "  C12 


+0  =  P0°  +  P2°  008  2cp 

where 


c  _  .  uo 
C20  "  c12  =  Y 


with 
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The  constraint  stresses,  which  will  occur  in  the  sheet  with  the  reinforced^ 
cutout,  are  included,  by  merely  adding  to  the  expressions  for  Pq°  and  Pg 

the  additional  terms  appearing  in  expressions  (3h) •  Thus, 


where 


in  which 


=  P0  +  P2  008  2<p 


c30loge  r 


~\ 

> 


+  C^2 


(21) 


The  constants  c^Q  and  c ■  in  expressions  (3h)  can  and  must  he  set 

equal  to  zero,  since  otherwise  they  would  give  constraint  stresses  which 

do  not  diminish  with  increasing  distance  r  from  the  cutout.  The 

constant  c^q  has  teen  omitted  since  it  would  not  contribute  anything  to  the 


stresses  or  displacements.  In  expressions  (21),  C£q  and  cj_2  are  given 
by  the  original  state  of  stress,  while 


?) 


'30' 


'22' 


and 


'42 


^equal  to 

are  the  unknowns.  These  unknowns  are  determined  for  a  given  reinforcement 
ring  by  putting  expressions  (21)  into  equations  (7a)  and  (7b),  which  express 


and  Ir,  and  by  satisfying  equations 
metric  term.  With  the  abbreviations: 


Ar 


(7a)  and  (7b)  in  every  trigono- 
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c30al  k30 


c22al  _  k22 


5  V 


\  El6(1  *  V>  ■  h 


A r  Er 


5  it'1  ♦  v)  - 


the  following  result  is  obtained  (see  appendix  D) 


0^(1  - 


1  21 
5 1  +  8  +  — 


30  1  +  OqI  2(1  +  V) 


1/  ,  .  ^  Av  -  16& 

tl1! ' 1  9aiU  +  v  . 


tv  -  16 


C1  +  35  + 


-  1  +  8  + 


1  +  V 


)  1  0^(1.+  48) 

fcc-^(6  -  88)  +  2] 


k42  =  ao 


|(l  +  36  +  ij  [1  -  0^(1  +  48)]  -  -1  +  8)  [ax(6  -  88)  +  2] 

C1  +  35  +  tK hrTr)  -  i1  + 8  +  rh)  [ai(6  - 86)  +  2] 
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Putting  expressions  (21)  into  equations  (4)  and  omitting  the  terms  in  c2q 
and  the  expressions  for  the  constraint  stresses  are  seen  to  be: 


Tr<pL  =  -  [&22(^)'4  +  2k^(^)“2  ] sln  2? 

/ 

In  order  to  calculate  the  stresses  in  a  plane  sheet  with  a  reinforced 
cutout  under  an  external  edge  loading  causing  an  original  state  of  stress 
(i.e.,  state  of  stress  if  there  were  no  cutout  at  all  in  the  sheet) 
given  by  ax  =  ao  =  Constant,  it  is  necessary  to  calculate  6 '  *1' 

and  <x^  in  accordance  with  expressions  (22)  from  the  dimensions  of  the 
reinforcement  ring,  determine  ^30*  ^22 ’  8X1(1  ^42  directly  from  equa¬ 

tions  (23),  find  the  constraint -stress  distribution  from  equations  (24), 
and  add  this  to  the  original  stress  distribution,  equations  (19)* 

The  case  of  uniform  tension  corresponds  approximately  to  the  experi¬ 
ments  carried  out  for  this  report  and,  with  the  analytical  calculation 
of  the  stresses  in  the  test  specimen  (appendix  G) ,  serves  as  an  illustrative 
example . 

Pure -bending  stress.-  If  an  I-beam  consisting  of  a  rectangular  web 
and  flanges,  as  shown  in  figure  1,  is  subjected  to  a  pure-bending  moment  M 
about  an  axis  z  perpendicular  to  the  plane  of  the  plate,  then  the  stress 
distribution  in  the  plate  will  be: 


M 


where 
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In  polar  coordinates  this  state  of  stress  is  given  "by: 


: — r(cos  39  “  cos  9) 

4lw 


\ 


M  \ 

- - r(3  cos  9  +  cos  39) 

%  4^ 


> 


Trcp  =  ~£f^r(sin  V  +  8in  3cp^ 


(25) 


It  can  he  verified  by  equations  (2)  that  the  Airy  function  for  this 
original  state  of  stress  is  given  by 


^o 


cos  cp  + 


COS  39 


where 


w 


By  substitution  into  equations  (7a)  and  (7b)  it  is  found  (see  appendix  E) 
that  it  is  not  possible  to  design  a  reinforcement  ring  producing  exactly 
zero  constraint  stresses. 

To  determine  the  approximate  additional  stresses  (in  this  case 
unavoidable)  due  to  a  cutout  reinforced  by  a  ring  of  given  dimensions, 
there  must,  as  in  the  preceding  case,  be  added  to  the  Airy  function  the 
additional  terms  in  expressions  (3b)  corresponding  to  the  same  trigo¬ 
nometric  orders  as  occur  in  Thus,  for  the  stresses  in  the  sheet 

under  pure  bending  with  a  reinforced  cutout 

f  cos  9  +  cos  3qp 

P1  =  c21r3  +  c31r_1 
P3  =  Cj_2r3  +  c23r”3  +  °43r 


(26) 
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The  constants  and  c-^  in  expressions  (3h).  have  heen  omitted  since 

they  would  otherwise  violate  the  assumed  condition  of  constraint  stresses 
decreasing  to  small  values  at  the  edges  of  the  sheet-  The  constant  c-^ 

has  heen  omitted  since  it  has  no  effect  on  the  calculations- 


Substituting  equations  (26)  into  equations  (7a)  and  (7b),  determining 
the  constants  c^,  c2^’  c43  (c2i  c13 >  aa  seated,  are  fixed 

by  the  applied  bending  moment)  so  that  equations  (7a)  and  (Tb)  are 
identically  satisfied  (see  appendix  F) ,  and  using  the  abbreviations 

c3  .  1 

ajt  k3 


c23al  k23 


c43al’3  S  k^3 
C3iai'3  =  k31 


Man 


Er 

Es 


+ 


v)  s  a 


Er 

Es 


—  24(1  +  V) 
a3t 


(27) 


and  the  assumption  V  =  0*3>  the  following  results  are  obtained: 
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With  the  constants  k^/  k^,  ^31  determined  ^rom  equations  (2.7) 

and  (28) ,  the  stresses  in  the  sheet  can  he  readily  calculated  by  adding 
the  constraint  stresses,  as  given  hy  equations  (29),  to  the  original 
stresses,  as  given  hy  equations  (25) . 

In  many  actual  cases,  the  reinforcement  rings  are  quite  narrow,  so 
that  if  the  approximations  6  =  0  and  i  =  0  are  made  only  negligible 
errors  will  he  introduced.  The  stresses  in  the  sheet  are  then  functions 
of  only  the  cross-sectional  area  of  the  ring.  This  is  true,  of  course, 
for  any  type  of  external  edge  loading  on  the  sheet  as  well  as  for  a  pure 
tending  load.  The  following  table  gives  the  values  of  the  constants  k 
for  pure  bending  for  different  values  of  the  area  concentrated  in  such  a 
line  reinforcement. 


PURE  BENDING 
[5  =  O',  i  =  0] 


> 

4* 

H 

OJ 

HI 

d 

0 

1/2 

1 

2 

VK 

3*00 

1.857 

1.48 

1.182 

Wr 

0 

0 

OJ 

1 

1 

00 

VJl 

-.486 

-.182 

2k3o/K 

-6. 00 

-2.31 

.426 

3.22 

Figure  2  shows  the  stress  distribution  in  the  sheet  for  a  =  1  and, 
for  comparison,  for  a  =  0  (i.e.,  no  reinforcement  around  the  cutout). 


Summary  of  Special  Cases 

Table  I  gives  a  brief  summary  of  the  numerical  results  of  the 
theoretical  investigation  of  the  special  cases  treated  here.  In  this 
table,  W.R.  denotes  the  ratio  of  the  weight  of  material  added  (in  order 
to  form  the  reinforcement  ring)  to  the  weight  of  material  cut  out-  Assuming 
that  material  of  the  same  density  is  used  for  both  ring  and  sheet,  this 
weight  ratio  is  given  by 


Ar  (3 
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Table  I  shows  also  the  ratios,  at  the  cutout 


of  maximum 


total  (i.e.,  original  plus  constraint)  stresses  to  maximum  original  stresses. 
In  the  column  headed  Remarks,”  the  term  "satisfactory”  means  that  the 
cross  section  concerned  has  practical  dimensions  and  that  at  the  same 
time  it  prevents  high  stresses  in  the  sheet. 


An  examination  of  table  I  shows  that  it  is  theoretically  possible  to 
obtain,  for  the  cases  of  hydrostatic  stress,  pure  shear,  single  tension 
(or  compression),  and  pure  bending,  zero  or  minimum  constraints  and  minimum 
total  stresses  by  reinforcement  rings.  However,  except  for  hydrostatic 
stress,  the  required  cross  sections  of  such  rings  would  have  very  high 
moments  of  inertia  but  very  small  widths  and  areas.  Such  sections  are 
extremely  difficult  to  design  in  practice.  The  table  also  shows,  on  the 
other  hand,  that  in  all  cases  considered  here,  not  the  minimum  but  at 
least  fairly  low  total  stresses  in  the  sheet  can  be  obtained  by  using 
appropriate  reinforcement  rings  of  practical  dimensions.  It  appears, 
moreover,  that  in  several  instances,  such  rings  may  weigh  less  than  the 
material  removed  from  the  sheet  to  form  the  cutout. 


TEST  ON  PLANE  SHEET  WITH  REINFORCED  CUTOUT 


For  the  experimental  part  of  this  research  a  test  was  made  on  a 
plane  sheet  with  a  reinforced  circular  cutout  under  a  tensile  load, 
realized  by  a  heavy  I-beam  transmitting  four  concentrated  loads  to  the 
sheet  (see  figs.  3  and  4)  «  The  loads  were  produced  by  adjustable  jack- 
screws  and  measured  by  calibrated  strain  gages  on  the  eight  connecting 
links. 

Before  carrying  out  the  main  test,  a  preliminary  test  was  made  with 
a  sheet  of  the  same  dimensions  but  without  cutout.  The  purpose  of  this 
preliminary  test  was  to  ascertain  whether  the  means  of  application  of 
the  load  would  produce  a  uniform  tension  stress  =  aQJ  *  Two  different 

total-load  stages  were  used,  with  the  loads  distributed  as  evenly  as 
possible  over  the  length  of  the  I-beam.  The  total  loads  for  the  preliminary 
test  were  1475  and  3016  pounds.  For  the  main  test  (i.e.,  with  cutout), 
the  loads  were  3000  and  4020  pounds. 

The  locations  of  the  strain  gages  (including  rosettes)  together  with 
the  results  of  experimental  measurements  in  the  sheet  without  cutout  are 
shown  in  figures  5  and  6. 

From  the  measured  values  of  the  test  on  the  sheet  without  cutout, 
it  is  seen  that  the  axial  strain  ex  was  distributed  fairly  uniformly 
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throughout  the  sheet-  At  the  higher  load,  the  maximum  deviations  from  the 
average  measured  value  were  -10  and  6- 4  percent  while ,  from  the  expected 

value  of  6  =  =  - 32^2. - -  =  2.98  X  10  the  maximum 

x  2btEs  o.052  x  18.5  X  10.5  X  106 

deviation  was  10.4  percent. 

Near  .the  loaded  edges  rosette  strain  gages  were  attached  to  show  the 
influence  of  the  lateral  constraint  of  the  riveted  Joint  of  sheet  and 
I-beam.  The  small  average  magnitude  of  the  transverse  strain  in 

the  sheet  near  the  loaded  edges  indicated  that  in  fact  the  load  was 
transmitted  by  the  heavy  I-bars  in  such  a  way  that  the  transverse  strain  ey 

^and  not  the  stress  was  practically  zero  at  these  edges.  If  the 

stress  cv  had  been  zero  at  the  edges,  then  the  magnitude  of  g _  would 

*7  J 

have  ‘been  much  higher,  namely  e  =  -Ve  =  -0*36  • 

1/  Jv  i 

Having  checked  this  type  of  original  stress  (or  strain)  distribution 
in  the  sheet  obtained  with  the  particular  loading  used  here  (see  fig-  3) > 
a  circular  hole  was  made  in  the  center  of*  the  same  sheet  specimen  and  was 
reinforced  ty  a  ring  of  dimensions  given  in  figure  3*  Strain  rosettes  and 
simple  strain  gages  were  then  placed  on  the  sheet  in  the  symmetric  posi¬ 
tions  indicated  in  figure  7*  Two  different  loads  were  applied,  3000  and 
4020  pounds.  Since  hoth  loads  produced  proportional  results,  only  those 
corresponding  to  the  higher  load  are  given  here.  The  results  of  the 
measurements  are  shown  in  figure  7* 

It  may  he  of  interest  to  compare  the  strains  measured  in  the  actual 
test  specimen  with  the  strains  calculated,  in  accordance  with  the  theory, 
on  the  basis  of  an  infinite  sheet.  Since  the  boundary  conditions  of  the 
test  piece  are  obviously  different  from  those  of  an  infinite  sheet,  the 
strains  may  be  expected  to  be  different  in  these  two  cases-  A  comparison 
of  the  strains  for  these  cases  may  nevertheless  be  instructive  in  indicating 
the  influence  which  the  finite-edge  conditions  have  on  the  strain  (or 
stress)  distribution  in  the  sheet.  This  comparison  is  given  by  the  strain 
diagram  in  figure  7.  Because  of  the  symmetry  about  both  x-  and  y-axes, 
the  theoretical  values  for  the  infinite  sheet  are  given  only  at  the  points 
indicated  in  the  diagram. 

From  figure  7,  it  will  be  observed  that,  qualitatively,  the  finite 
sheet  behaves  quite  similarly  to  the  infinite  sheet.  For  example,  in 
both  cases  the  transverse  normal  strain  f€  \  decreases  steadily 

'  9'— =1.38 

ai  ■ 
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with  cp  from  a  maximum  at  cp  =  0°  to  a  minimum  at  cp  =  90°. 

Also  (e •  \  near  the  rigid  bars  varies  in  "both  cases  from  a 

Vx^x=9-2  in. 

maximum  at  the  edge  of  the  sheet  to  a  minimum  in  the  center*  (See  graph 
"below  strain  diagram,  fig*  7*)  The  radial  strains  were  all  of  a  relatively 
low  order  of  magnitude,  "but  it  can  "be  stated  that  here,  also,  "both  the 
order  of  magnitude  and  the  variation  of  €  corresponding  to  the  calcula¬ 
tions  for  an  infinite  sheet  were  similar  to  those  of  the  test  specimen. 


The  effect  of  the  particular  "boundary  conditions  of  the  test  piece 
appears  to  he  pronounced  at  two  places.  The  maximum  transverse  normal 
strain  ^  (perpendicular  to  the  radius  and,  at  cp  =  0,  in  the  direction 

of  the  general  tension  stress)  near  the  cutout  ring  (  —  =  1. 3© )  is  higher 

Vai  / 

for  the  test  piece  than  for  the  infinite  sheet,  the  percentage  difference 


being 


5.06 


-  4.07 


X  100  =  19*6  percent . 


The  axial  strains  near  the  loaded 


5.06 

bars  (see  graph  below  strain  diagram,  fig.  7)  decrease  more  sharply  in  the 
actual  specimen  than  in  the  specimen  calculated  as  an  infinite  sheet,  the 


local  percentage  difference  at  the  center  being 


3.17  -  2.32 
2.32 


X  100  =  36*  6  per¬ 


cent,  although  the  analytical  values  were  practically  equal  to  the  experi¬ 
mental  values  at  the  quarter  points  of  the  sheet*  The  difference  at  the 
center  must  be  due  to  the  fact  that  in  the  actual  test  piece  the  rigid 
bars  caused  constant  axial  displacements  and  zero  transverse  displacements 
at  the  finite  loaded  edges,  whereas  in  the  theoretical  work  an  infinite 
sheet  was  treated  with  constant  axial  stress  at  the  remote  loaded  edges* 


The  maximum  strain  measured  in  the  sheet  was  the  shearing  strain 

at  cp  =  45°,  and  it  will  be  observed  that  this  value  (7*66  X  10  was 
practically  unaffected  by  the  difference  in  the  boundary  conditions  between 
the  actual  and  the  theoretical  specimen,  since  the  experimental  and  the 
theoretical  values  are  practically  equal. 


The  calculations  for  the  infinite -sheet  specimen,  based  on  the  formulas 
developed  in  the  analysis  preceding  the  experiments,  are  shown  in  Appendix  G* 


It  may  be  remarked  that  in -these  calculations  account  was  taken  of 
the  rivet  holes  in  the  cutout  ring.  These  holes  had  the  effect  of  reducing 
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the  actual  cross-sectional  area  of  the  ring.  To  determine  the  amount  of 
this  reduction,  a  test  with  the  same  size  and  spacing  of  the  rivet  holes 
as  in  the  sheet  specimen  was  carried  out. 

The  ratio  of  the  effective  area  to  the  full  geometric  area  was  thus 
found  to  he  0.80.  Therefore  in  the  calculation  of  the  theoretical  stresses 
in  the  test  specimen  used  in  the  experiment  the  cross-sectional  area  of  the 
cutout  ring  was  taken  as  0.80  times  the  full  geometric  area. 


CONCLUDING  DISCUSSION 


The  problem  first  considered  was  what  must  be  the  cross  section  of 
a  ring  reinforcing  a  circular  cutout  in  a  plane  sheet  in.  order  that  the 
stresses  in  the  sheet  remain  unchanged  by  the  cutout.  The  general  solution 
to  this  problem  is  given  by  equations  (7a)  and  (7b)  •  In  these  equations 
the  required  distributions  of  moment  of  inertia  Ir  and  of  cross-section 

area  Ap  along  the  circumference  of  the  reinforcement  ring  are  expressed 
in  terms  of  the  stress  function  (with  coefficients  Fn  and  Qn  as 

defined  by  equations  (3a)  and  (3b))  for  the  original  stress  distribution, 
that  is,  for  the  stresses  in  the  sheet  without  cutout  or  reinforcement 
rings. 

It  was  found  that  these  expressions  in  terms  of  the  circumferential 
angle  cp  for  Ir  and  Aj.  lead  to  physically  possible  (i*e-,  positive) 

values  only  in  a  limited  number  of  cases  (for  example,  when  the  original 
stress  function  has  no  or  only  one  trigonometric  term,  as  in  centric 
symmetry  and  in  pure  shear,  respectively) •  In  the  other  cases,  which 
include  uniaxial  tension  and  pure  bending,  constraint  stresses,  that  is, 
additional  stresses  due  to  the  reinforced  cutout,  are  unavoidable.  A 
method  of  calculating  such  stresses  in  the  sheet  for  a  given  (constant) 
cross  section  of  the  reinforcement  ring  was  developed,  based  on  the 
requirement  that  the  constraint  stresses  and  displacements  vanish  suffi¬ 
ciently  rapidly  with  increasing  distance  from  the  cutout  so  as  to  have  a 
negligible  influence  at  the  edges  of  the  sheet.  The  formulas,  which  are 
straightforward  and  convenient  to  apply,  were  derived  in  detail  for  the 
cases  of  pure  shear,  uniform  axial  tension  (or  compression),  and  pure 
bending.  For  example,  in  the  case  of  uniform  axial  tension,  it  is 
necessary  merely  to  calculate  the  values  of  the  dimensionless  constants  ip 
and  ot^  (proportional,  respectively,  to  the  cross-sectional  moment  of 

inertia  and  the  area  of  the  reinforcement  ring)  from  the  given  data  in 
accordance  with  the  definitions,  equations  (22),  of  ip  and  o-^.  The 

vaPues  of  the  constants  k^Q,  k22>  and.  k^g  follow  from  the  elastic 

properties  of  the  ring  and  are  determined  by  equations  (23) .  The 
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constraint  stresses  (which  must  be  added,  to  the  original  stress  distri¬ 
bution  Oy.  =  aQ  =  Constant)  then  follow  readily  from  equations  (24). 

In  the  formulas  derived  for  the  cases  deemed  technically  important, 
it  can  be  seen  that  for  the  narrow  reinforcement  rings  commonly  used  only 
smifli 1  errors  will  be  introduced  in  the  numerical  calculations  if  the  moment 
of  inertia  is  taken  as  zero*1  The  results  of  the  analysis  therefore  show 
that  the  stresses  in  a  plane  sheet  with  a  cutout  are  a  function  only  of 
the  cross-sectional  area  of  the  ring  reinforcing  the  cutout,  while  the 
moment  of  inertia  is  of  practically  no  influence.  This  means  that 
the  ring  experiences  primarily  tensile  or  compressive  stress  result¬ 
ants  and  not  bending  moments. 

In  general  the  approximation  given  by  the  formulas  of  the  report  will 
be  closer  the  smaller  the  constraint  stresses  are  at  the  edges  of  the  finite 
sheet . 

In  the  experimental  part,  a  test  was  made  of  a  plane  sheet  with  a 
reinforced  cutout  subjected  to  a  uniform  tensile  displacement.  The 
purpose  of  this  experiment  was  to  see  how  great  an  error  is  produced  by 
boundary  conditions  which  differ  from  those  for  which  the  theoretical 
formulas  are  exact-  In  particular,  the  sheet  specimen  was  (of  course) 
not  infinite,  while  at  the  loaded  edges,  the  axial  displacements,  and  not 
the  stresses,  were  constant  (with  zero  transverse  displacements  there) • 

On  the  two  other  opposite  edges  the  normal  and  the  tangential  edge  tractions 
were  zero. 

It  was  found  that  qualitatively  the  strain  (and  therefore  stress) 
distribution  in  the  sheet  specimen  was  quite  similar  to  that  predicted 
by  the  theoretical  formulas  for  an  infinite  specimen  with  the  same^rein- 
forced  cutout.  Quantitatively,  the  chief  effect  of  the  actual  boundary 
conditions  seemed  to  be  at  the  loaded  edges  where  the  axial  strains 
decreased  more  sharply  toward  the  center  of  the  sheet  than  predicted  by 
the  formulas  for  the  infinite  sheet- 


The  transverse  strain 


’qp 


at  the  transverse  center  line  of  the 


sheet  near  the  cutout 


was  actually  ab'out  20  percent 


^This,  of  course,  is  true  only  if  there  is  no  concentration  of  the 
original  stress  along  any  radius  of  the  cutout  circle,  a  case  which  appears 
for  a  concentrated  load  at  an  external  edge  but  which  was  not  considered 
in  this  report. 


32 


NACA  TN  No.  1852 


higher  than  according  to  the  theory  for  the  infinite  sheet.  In  other 
respects,  however,  the  quantitative  results  appeared  not  to  he  greatly 
affected  by  the  particular  boundary  conditions  of  the  test  specimen, 
s^nce  the  experimental  values  of  the  strains  were  fairly  similar  to  the 
analytical  values  based  on  an  infinite  sheet. 


Polytechnic  Institute  of  Brooklyn 

Brooklyn,  N.  Y-,  June  27,  1947 
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APPENDIX  A 

STRESS -STRAIN  RELATIONS  OF  REINFORCEMENT  RING 


The  theory  of  curved  teams  states  that  the  change  of  curvature  Kr 
of  a  ring  due  to  a  tending  moment  G  is 


G 

E  I 
r  e 


where  the  effective  moment  of  inertia  I0  is  given  ty 


and  where  r  =  rQ  denotes  the  radius  of  the  neutral  fiber,  while  t  is 

the  width  perpendicular  to  the  plane  of  bending  of  a  cross  section,  and  r0, 
moreover,  is  given  ty  (cf.,  for  example ,  reference  l) 


The  following  table,  however,  shows  that  for  a  cross  section  of  I-shape 
(for  example)  with  values  of  d/a  up  to  0-4,  a  ring  may  be  treated  as  a 
straight  beam,  so  that  it  is  permissible  to  put 


d 


I  =  I 
e  r 
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COMPARISON  BETWEEN  CURVED  BEAM  AND  STRAIGHT  BEAM 


d/a 

d 

al  “  2 
( in. ) 

ro 

(in. ) 

(in.1*) 

(in.4) 

0.2 

6 

5 - 973049 

0.041855 

0.04187 

•3 

6 

5-93044 

.119638 

.124966 

.4 

6 

5-865331 

.252761 

.260266 

It  follows  then  that 


Kr  = 


E  I 
r  r 


as  In  equation  (6b).  The  extension  strain  er  and  the  curvature  **r  of 
a  ring  are  given  in  terms  of  the  radial  and  tangential  displacements  urQ 


and 


"t 


at  the  neutral  fiber  r  =  rD  by  the  following  well-known  relations: 


6  a  =  u-i-  +  Uv, 
r  ro  ro 


Vs  - '(% + 

Since  the  radial  displacement  Up  of  the  ring  for  any  given  cp  will 
be  constant  along  the  width  of  the  ring,  it  follows  that 

*ra2  =  +  Up)  =  +  Up) 

y  r=r  \  x  x/r=a1 
o  1 

The  tangential  strain,  however,  will  vary  along  the  width  as,  in  fact, 
can  he  seen  from  figure  8, 


a- 


..  d 
"r  2a 
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i 


Hence 


Thus,  equations  (6a)  and  (6h)  follow. 
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APPENDIX  B 


EXPLICIT  EXPRESSIONS  FOE  STRESS  EESULTANTS  AND 


STRAINS  OF  EEINFOBCEMENT  EING 


From  figure  1  the  following  equilibrium  conditions  between  the  stress 
resultants  T,  G,  and  N  in  a  ring  and  the  unit  load  forces  Y  and  Z 
can  be  derived: 


T  +  N  =  -Ya1  =  - 


NL 


N  -  T  =  -Za  =  -(cr)  a-,t 
V  /  r=a. 


G  -  NrQ  =  -Ya^e  =  - 


C  r(p)r=i 


Noting  that  rQ  £  a,  e  SJ  —  (see  appendix  A),  the  last  equation  can  be 


written  as 


G  -  Na  =  -T^t 


Equations  (Bl)  to  (B3)  can  “be  solved  for  T  and  G  as  follows 
From  equations  (Bl)  and  (B2)  , 


T  +  T  =  a 


lt(°r  “  Tr<p)  r=« 


Moreover,  from  (Bl)  and  (B3) > 


G  =  ~T&  a/t  f  (Tr„)  d„ 
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By  putting  ar  and  t  in  terms  of  the  stress  function  (relations  (2) ) 
equation  (B4)  can  "be  directly  integrated  with  the  results: 

for  n  5/  1, 


T  =  a^  (r-1V)r=ai  +  c2  ain  9+03  cos  cp 


G  =  alt  §  (r  1^)r=a1"  al2t  (r_2llr)r=a1  +  C1  “  c2a  sln  "  °3a  008  <P 


for  n  =  1, 


T1  ~  c2 


sin  cp  +  c 


3 


cos  q> 


G1  “ 


c2a 


sin  cp 


-c3a 


COS 


<P  +  C-JL  +  ax2t  (r-1V  •  -  r“2ilf)r=ai 


Substitution  for  t  by  means  of  relations  (3)  will  give  expressions 
for  T  and  G  in  terms  of  Pn  and  cp. 

Expressions  in  terms  of  Pn  for  the  extension  strain  and  curvature 

can  be  obtained  by  putting  equations  (2)  into  equations  (5)  and  integrating 
the  latter  to  find  ur  and  u-fc.2  Thus, 


2 

Two  arbitrary  functions  fj_(r)  and.  f2(cp)  will  appear  as  a  result 

of  this  integration.  However,  it  can  be  shown  that,  for  compatibility 
between  deformation  values  (involving  also  the  shear  strain  7rc p)  and 

the  equilibrium  conditions,  the  functions  must  have  forms  which  are 
already  included  in  the  general  expressions  (3a)  and  (3b)  for  the  Airy 
function. 
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Es  (^r  +  *r)  = 


f'+  2 V  +  t 
„2 


dr  +  p  (t  +  10  -  v(+  +  t) 


E 


s 


By  using  the  expression  (3a)  for  ty,  the  strains  in  terms  of  P. 
are  obtained. 


and  cp 
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APPENDIX  C 


IMPOSSIBILITY  OF  ZERO  CONSTRAINTS  IN  CASE  OF  UNIFORM  TENSION 


=  Constant^ 


From  the  expression  (20)  for  i|r  , 


& 


Op 


po  =  fr 


Po  -  -2r2 
r2  ~  k 


Substituting  these  values  into  (7a)  and  (7b),  putting  c2  =  =  0  (for 

symmetry  about  both  x-  and  y-axes),  and  simplifying  give  the  following 
expressions  for  the  distribution  of  cross-sectional  area  and  moment  of 
inertia: 

Ea 


Ar  =  at 


1  +  cos  2<p 


Er  (1  -  v)  +  |l  +  V  +  kb(l  +  v)]  cos  2qp 


-  Es  (1  -  &)  (1  +  cos  2cp) 

Iv.  =  a3t  m2 


ci 

a^a 


2(1  -  v)  +  6(1  +  v)  cos  2cp 


where  6  =  d/2a- 


It  is  not  difficult  to  see  that  Ir  can  he  made  positive  for  all 

values  of  cp  (h y  choosing  a  proper  value  for  c^)  hut  that  neces¬ 

sarily  becomes  negative  and  Infinite  for  some  values  of  cp*  Hence  a 
perfectly  equivalent  reinforcement  is  not  realizable  In  this  case. 
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V. 


2  +  k30  al 


gp(i  -v) 

2(1  +  V)  30 


D 


2  “ 


+  k22(6  -  85)  +  k^2  ^ 


K 


•  _  r  1  -  v  k3o 

l  =  1i  [12(1  +  v)°o  "  6  _ 


°o.-  .  /.  -_v  .  .  ro 


“(1-8)  +  k22(l  +  35)  +  k42(l  +  5)  =  i±(^  -  k22 


4v  ~  165\ 
1  +  v  / 


r 


The  first  of  these  four  equations  can  he  solved  immediately  for  k^Q*  The 

third  equation  determines  the  value  of  Kj_,  which,  however,  has  no  influence 

on  the  stresses  in  the  sheet.  Hence  this  equation  need  not  he  considered 
any  further  here*  The  second  and  fourth  equations  can  he  solved  simul¬ 
taneously  for  k22  and  T^e  resu^-^s  are  *>y  equations  (23). 
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APPENDIX  E 

PRACTICAL  IMPOSSIBILITY  OE  PERFECT  EQUIVALENCE  IN 
CASE  OF  PURE  BENDING 

From  the  expression  for  it  is  seen  that,  corresponding  to  the 

(original)  state  of  stress  due  to  pure  bending, 

P1  =  C21r3 


P3  '  °13r 


<=21  =  3c13  -  -8ll 


Substitution  into  equations  (7a)  and  (7b)  leads  to  the  following  expressions 
for  Ap  and  Ir  =  c2  =  0,  "because  of  the  antisymmetry  about  the  x-axis 

and  the  symmetry  about  the  y-axis): 


cos  cp  +  cos  3cp 


at _ 3aj  tc13 

(2  +  95) (1  +  v)  6  -  2V  +  5(3V  -  1) 

(2  +  95) (1  +  V) 


cos  qp  +  cos  39 


ir  -  -s2s-u  -  25)  Hs^sa 

24(1  +  V) 


6(1  +  6)  /  (1  -  26)  f  cos  cp  +  cos  3cp 


cos  3<P 


It  will  be  shewn  now  that  Ay  and  Ir  for  zero  constraints  can  both  be 

made  positive  for  any  cp  only  if  the  value  of  the  ratio  6  (identical 
with  d/2a)  is  extremely  low,  too  low  in  fact  for  practical  purposes. 
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^3 


4 


% 


From  the  second  of  the  two  equations,  in  order  that  Ir  he  everywhere 
positive  it  is  necessary  that 

°3  -  6(1  +  8) 

“1  t013  ' 

Substituting,  then,  for  into  the  expression  for  the  area  Ar  gives 

a-k  2(l  +  8)  cos  qj  +  cos  l<p 

Ar 


(2  +  98) (1  +  v) 


6  -  2V  +  8(^V  -  1) 
(1  +  v)(2  +  98) 


f 


cos  qp  +  cos  39 


The  value  of  this  expression  for  Ap  will  he  positive  regardless  of 
and  only  if 


2(1  +  5) 


6  -  2v  +  8(3V  -  1) 
(1  +  v)(2  +  95) 


cp  if 


With  V  =  0»3>  the  positive  root  of  this  quadratic  equation  in  5  is 


5  =  0.006814- 


In  practical  design,  this  low  value  for  5  (identical  with  d/2a,  fig.  l) 
will  he  incompatible  with  the  large  required  moment  of  ' inertia 

a3t 


Ir  24(1  +  v) 


(1  -  28) 


DERIVATION  OF  EQUATIONS  (28)  PERTAINING  TO  CASE  OF  PURE  BENDING 
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The  second  and  fourth  equations  are  linear  equations  in  k^  and  k^, 

which  can  therefore  he  readily  solved  for  these  two  unknowns.  By 
eliminating  k^  from  the  first  and  third  equations,  k^  can  be  solved 

for  in  terms  of  kg^  and  k^.  The  results  are  given  by  equations  (28) 
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APPENDIX  G 

CALCULATION  OF  STRAINS  IN  INFINITE  -  SHEET  SPECIMEN  WITH  SAME 
REINFORCED  CUTOUT  AS  ACTUAL  SPECIMEN 


The  dimensions  of  the  specimen  are  given  in  figure  3 
dimensions, 

a  =  3*66  inches 


From  the 


a-,  =3*5  +  -^-  =  3- 8l  inches 
^  16 


t  =  0-052  inch 

1. 

5  =  —  =  - — -  =  0-0426 

2a  2  X  3-66 


Ae  =  0.80  X  0.434)  =  O.IO85  square  inch 

Ae(1  +  =  0.1085  X  1.3  _  0>7l^ 

at  3.66  X  0.052 


^(j) 


X  0.434  .3  u 

-  =  1.1  x  10  inch4' 


12 


ix  ^  6(1  +  V)  — 


7-8  x  1:1  X  10-3  -  3-37  x  10-3 


a3t  (3 .  66r  x  0.052 


From  equations  (23)  substitution  of  the  preceding  values  leads  to  the 
following  numerical  results  for  the  constants  k22,  and  k^: 
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k30  "  _0*173cro 


k22  =  °-022®ao 


\2  ~  •°'252oo 


From  equations  (24)  it  then  follows  that  the  calculated  constraint 
stresses  in  the  sheet  will  he: 


*1 


/  VI 

\“2 

_  /~r>  \  “4 

\  -2~T 

-0,173(f 

0  ■ 

-  1.008^  1 

^  +  0.1368^^  ^  cos  2cp 

T 

_ 

rq>,_ 
co  ” 

0.1368 

©“  -  “■»*©■* 

sin  2qp 

cos  2qp 


The  calculated  stresses  without  cutout  (see  equations  (19))  are: 

a0 

a  -  — (i  -  cos  2cp) 
ro  2 


0<po  =  ~2^1  +  cos  2cP) 


Trcpo  =  T  sln  2(P 


Hence,  adding  the  constraint  to  the  original  stresses,  the  total  stresses 
in  the  sheet  at  any  point  (r,  cp)  are  found  to  "be: 


Gr  =  ao  < 


0-500  -  0.173 


acp  =  ao  '  0,500  +  0,173 


t  =  a 
rep  o 


(qV2  -  o.  500  +  0.1368  (i)'4  -  i-°o8^)'2J 

1368(f1)-‘ 

/rV2  / 

0.500  +  o-5o4(_-J  .  0.1368b 


cos  2cp 


0- 500  +  0* 


r  \ 

Fx 


cos  2cp 
sin  2cp 


h8 
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From  these  stresses,  the  strains  expressed  in  polar  coordinates  at  any 
point  are: 


2(1  ,t_V)T 

'rq>  Eg  rep 


The  strains  ex  in  the  direction  of  the  load  were  obtained  by  first 
determining  the  Cartesian  stresses  o  and  a  according  to  the  relations 


1 

°x  2 


ay  =  2 


^r  +  "  (0r  _  a<p)  008  2(p  +  ^rep  sin  2cp] 

(ar  +  °<p)  +  (°r  *  acp)  008  2cp  '  2rrep  sin  2<PJ 


The  strain  ex  is  then  readily  obtained  by 

■  ij*  -  %) 

The  value  of  the  original  stress  aQ  without  cutout  is 


iff 


°°  2bt 


where  2b  is  the  width  of  the  sheet.  For  the  test  specimen  and  for  a 
load  of  h020  pounds. 


a  _  -  =  hi 60  pounds  per  square  inch 

0  18.5  X  0.052 
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For  the  strain  this  can  he  expressed  as: 


4l6o 


xo  10-5  X  106  x  4.02 


=  0-99  x  10 


-4 


The  actual  measurement  gave  1.04  x  10~\  that  is,  an  average  inaccuracy 
of  5  percent  in  the  strain  gages. 

The  numerical  results  of  the  preceding  calculations  are  given  in  the 
data  diagram  (fig.  7)  • 


50 


NACA  TN  No.  l8^>2 


1-  Timoshenko,  S. : 

and  Problems. 


EEZERENCE 


Strength  of  Materials.  Part  II  -  Advanced  Theory 
Second  ed. ,  D.  Yan  Nostrand  Co.,  Inc-,  19^-1*  ch.  II. 


NACA  TN  No.  1852 


51 


TABLE  I 

SUMMARY  OF  THEORETICAL  INVESTIGATION  OF  SPECIAL  CASES 


Loading 

A 

2a 

Jr 

at 

Jr. 

a^t 

W.R. 

(1) 

Ratios  of  maximum 
total  stresses  t6 
maximum  original 
stresses  at  r/a-^  =  1 

Remarks 

(3) 

°r 

°ro 

X 

Trq> 

Trcp0 

°*o 

(2) 

Hydrostatic 
or  =  Constant 

0 

1.1+3 

Any 

2.86 

1.00 

1.00 

1.00 

2 

Zero  constraints; 
satisfactory 

.1 

1-^3 

Any 

3.04 

1.00 

1.00 

1.00 

.2 

1-43 

Any 

3-21 

1.00 

1.00 

1.00 

Pure  shear 

t  =  Constant 
xy 

.2 

.1)28 

0.1025 

_ 1 

•  0875 

1.0 

1.0 

1.0 

4 

Zero  constraints,  tut  Ip 

too  high  to  conform  with 

Aj,  and  d/a 

.1 

.385 

.0032 

.457 

-.46 

1.30 

1.58 

Satisfactory 

.1 

.550 

•H55 

.865 

1.0 

1.0 

1.0 

Zero  constraints,  tut  Ir 

too  high  to  conform  with 

Ap  and  d/a 

Single  uniform 
tensile  (or 
compressive) 
stress 

ax  =  Constant 

.1 

.882 

•115 

1.69 

.822 

.872 

1.31 

Low  constraints,  tut  Ir 

too  high  to  conform  with 

Ar  and  d/ a 

.2 

•795 

.1025 

1.232 

.783 

•790 

1.43 

3 

Low  constraints,  tut  Ir 
too  high  to  conform  with 

A^  and  d/a 

0 

.615 

0 

I.23 

.750 

1.44 

1.63 

Satisfactory 

.1 

.461 

0 

.682 

.512 

1-31+ 

1.59 

Satisfactory 

Pure  tending 

0 

.385 

.0320 

•770 

1.015 

■995 

1.015 

2 

Almost  zero  constraints,  tut 
Ir  too  high  to  conform 

with  Ap  and  d/a 

0 

•385 

0 

•770 

1.01 

1.10 

1.34 

Satisfactory 

W-E*  ratio  of  the  weight  of  material  added  (to  form  reinforcement  ring)  to  the  weight 
of  material  cut  out. 

2 

No  reinforcement. 

3  ^  11 

Satisfactory  means  that  the  cross  section  has  practical  dimensions  and  that  at  the 
same  time  it  prevents  high  stresses  in  the  sheet.  "Zero  constraints"  means  that  the  rein- 
forcement  ring  is  exactly  equivalent  to  the  portion  of  the  sheet  which  has  teen  cut  out. 


polar -coordinate  system. 
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Figure  3.-  Experimental  setup. 
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Figure  4.-  Test  specimen, 
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369  364  360  369 

746  764  729  746 


363  379  372  374 

760  770  749  769 


Figure  5.-  Measured  strains  ex  x  10 ^  in  direction  of  loading.  Top 

values  are  for  a  total  load  4P  of  1475  pounds;  bottom  values  are  for 
a  load  of  3016  pounds.  Rosette,  A  ;  simple  strain  gage,  .  Drawing 
not  exactly  to  scale. 
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r4nrt“ir 

m5.0 1  3.77  2.32 
th  (4.26)  (3.84)  (3.17) 


I  4.  13 
32.6"  -.80 

!  .28 


5 .06  //  > 

m  //  /V 


7.75  — 
I  .59 

7.S 


.504(- 0.295) 
.61(1  .62 X 
.27  (OK/ 


.72 

.20 


\l ,94( I .88) 
,595(.935) 
K69  (7.62) 

iT06  \  (4.07) 
.094,  (.242) 


3.91  (4.00) 
-.79  (-.488) 
.25  (0) 


th 


(£x)x  =  9  2"  against  y 


Figure  7.-  Strain  measurements  on  plane  sheet  with  reinforced  cutout. 
Total  load,  4020  pounds;  rosette,  A  ;  simple  strain  gage,  p) ;  values 
in  parentheses  are  corresponding  theoretical  values  for  reinforced 
cutout  but  for  an  infinite  sheet  (see  appendix  G).  For  each  rosette, 
the  top,  middle,  and  bottom  values  are,  respectively,  normal  transverse 

normal  radial,  and  shear  strain  times  104.  For  each  simple  gage, 
values  given  are  those  of  axial  strains  e  v  x  10/ 
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